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Abstract. The vector radiative transfer equation is used for the solution of numerous problems 
in the field of random media optics. It describes the transformation of the Stokes vector of a light 
beam due to its propagation and scattering in a medium. However, the Stokes vector depends 
on the choice of a coordinate system. Thus, one needs to account for the rotation of coordinate 
systems within the framework of the standard formulation of the vector radiative transfer equation. 
The coordinate-free approach has advantages in the case of complex random media, including 
anisotropic and chiral ones. The aim of this paper is to derive the coordinate-invariant tensor 
radiative transfer equation. The derivation is based on the concept of the light beam tensor. We 
also study the relationship between the tensor and vector forms of the radiative transfer equation. 


1. Introduction 


The propagation of light in turbid media is often described within the framework of vector 
radiative transfer theory [1-3]. The vector radiative transfer equation [1—3] describes the 
change of the Stokes vector S of a light beam due to its propagation and scattering in a 
random medium. The components of the Stokes vector 7, Q, U and V can be related to the 
components of the electric vector E=E itt + Ezz of a plane monochromatic wave by the 
following equations: 


I = |E? + |E2|7 Q = |E? — | E2}? 


1 
U = 2Re (E, E3) V = —2 Im (E, E3) Z 


where the values of E; and E, are components of the electric field in the plane, perpendicular 
to the direction of propagation. We have omitted the common multiplier in equations (1) for 
the sake of simplicity. It should be pointed out that equations (1) only hold for monochromatic 
waves. In the more general case of non-monochromatic waves, the phases and amplitudes of 
the electric vector are not constant and the parameters (1) should be averaged with respect to 
time t. 

One can see that the Stokes vector is defined within the framework of the special coordinate 
system, which is attached to the direction of light propagation. This direction is changed many 
times during the process of light scattering in a random medium. This feature is accounted 
for by using special rotation matrices [2,3]. It is of general importance to derive the radiative 
transfer equation in a tensor form, which does not depend on the specific coordinate system. 
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The covariant methods were introduced into the optics of anisotropic and gyrotropic media 
by Fedorov [4, 5] and led to rapid progress in the field of anisotropic media optics [4—6]. This 
was due to the fact that the direct manipulation of vectors, dyadics and their invariants facilitates 
solutions, condenses exposition and provides results of greater generality. 

However, the coordinate-free approach has not been used much in the optics of random 
anisotropic media so far. The aim of this paper is to derive the covariant tensor radiative 
transfer equation and bring attention to covariant methods in radiative transfer theory. 


2. Light beam tensor 


Let us define the light beam as a superposition of different incoherent simple uniform waves. 
These waves can have random phases and different states of polarizations. However, they are 
characterized by the same speed and direction of propagation. The general case of incoherent 
light beams may be studied using the notion of the light beam tensor [5] 


Ra yore nE®O = iE =0 (2) 
AY 


where E = E |X+ E2y+E3Z is the electric vector, Se denotes the summation over all incoherent 
simple waves in a beam, and the vector 7 defines the direction of propagation. The dyadic 
notation F = E E* means 


EE} EE E\E% 
F = EE% EE% E2 E% . (3) 
EE} E3E% E3E% 


One obtains within the framework of the coordinate system, attached to the direction of 
propagation 7||z that E; = 0 and the tensor F (3) reduces to the density matrix [3,5] 


IE? Ey Es 
= : 4 
P=\ met [Ey 


The elements of the matrix p can be expressed in terms of components of the Stokes vector 
S (see equation (1)) 


ol I+Q0 U-iV i 
DON u+iv 1-ọ J (a) 


Thus, both (1) and (4) depend on the choice of the coordinate system. One can find the radiative 
transfer equation for the density matrix p in [3]. 
Fedorov [5] introduced the following invariants of the light beam tensor F: 


I=F, K = (F°), M =i(n*F), L= (FF*), (6) 


where t denotes the trace (e.g. F, = | E,|* +|E2|* +| E3|?, see equation (3)) and n* is the tensor 
with components 


ne. = CabcNb- (7) 


Here values of n, are components of the vector n and values eabe compose the Levi-Civita 
tensor [5] with components e123 = €23; = €312 = l and e213 = e132 = €321 = —1. Other 
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components of the Levi-Civita tensor are equal to zero. Equation (7) can be written in the 
following explicit form: 


0 —n3 N2 
n* = n3 (0) Ny 
=n? ny 0 


It follows from equation (6) that 
I=J+h+h 
K = J? + J2 + J +2 (FF + Fiz Fa + Fo3 F32) 
L= J? + J2 + J? + FÈ + FA + Fh + FA + Fat Fe 
M = i[n3 (Fi2 — Fo) +m (F31 — Fi3) +n, (Fos — F32)] 


where F;; = )), E; E3*, J; = Fj; and $`, denotes the sum over rays in a beam. 
It follows that 


I=])+h K = J? + J2 +2Fy Fz 
L= J? + J2 + F, + FA M = 21Im(F,2) 


at n||Z(n1 = n2 = 0, n3 = 1). One can see that the first invariant describes the intensity of 
the light beam J. The invariant M at ||Z coincides with the fourth component of the Stokes 
vector V (see equation (1)). Thus, the first and last components of the Stokes vectors (see 
equation (1)) are coordinate independent. It should be pointed out that this is not a case for 
components Q and U. 

The ratio K /I* determines the degree of polarization [5] 


2K 
p=\—-! (8) 


and the invariant L can be used to find the semi-axes of the polarization ellipse [5] 


a= J3;U+VL) b= J5U-VL). (9) 


For instance, it follows that 


a=b= 1/2 (10) 


at L = 0. This is the case for circularly polarized light [5]. 

Values (1) can be found from components of the tensor F. Thus, one can see that the 
light beam tensor and its invariants can be used for the complete description of the polarization 
characteristics of a light beam. 


3. The radiative transfer equation 


Let us now find the transfer equation for the light beam tensor in a random medium. It follows 
on general grounds that the change of the light beam tensor dF (7) in the direction, specified 
by the vector ñ, is due to two processes, namely, due to light scattering (dF) (n)) on the path 
dn from all directions to the direction 7, and light extinction (dF ) (#)) on the path dn due to 
beam propagation in a light scattering and absorbing medium 


dF = dF) +dF. (11) 
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Let us now calculate these contributions for a cylindrical volume with a unit cross section. 


First of all we note that light scattered by a single particle can be presented as a spherical wave 
> ., exp(ikR) 
E sca ( Pp ) © R 


where Ei (p) is the vector amplitude of the spherical wave scattered in the direction p, R is 
the distance to the observation point and k is the wavenumber. The amplitude E,ca(p) can be 
presented in the following form due to the linearity of the Maxwell equations: 


Esal P) = f, DEO) a2) 


where Eo (q) is the electric field of the incident wave, propagating in the direction specified 
by the vector g, and f (p, q) is the scattering tensor. Let us introduce the light beam tensor of 
the incident beam 


F= E @ES @ (13) 
and the light beam tensor of the scattered light 

F' =) EQS@BES P) (14) 
where the symbol s represents a single ray in a beam. It follows from equations (12)—(14) that 


FO) =) EOE SS PEE, fF ator (15) 
where f* = f * and f is the transpose ( Îr = fki) tensor. This formula establishes the law of 
transformation of the light beam tensor due to the scattering process from the direction g to 
the direction p. It is clear that the total light scattering in the direction specified by the vector 
n on the path dn from all directions n’ will differ from equation (15) due to the integration on 
the solid angle dQ;, namely 


dF (n) = Nan | dq FELDES ED 


where N is the number of scattering events in a unit volume. 
The removal of photons from the direction 7 on the path dn due to the extinction process 
can be described by the operator L 


dF (i) = —NLF() dn. (16) 
Thus, it follows from equation (11) that 


dF (n) F > > 3 >) +721 > > 
a= -NÊF@ +N | a2 fG n)F (a) f* R, A) +B) (17) 


where the term B(n) describes the internal emitting sources. This is the radiative transfer 
equation written within the framework of the coordinate-free approach. 

Let us now determine the operator L. The tensor F can be presented in the following 
form: 


Fm) = F.(n)5(n — no) + Fa) (18) 


where 6 (n — no) is the delta function. Values F.(n) and F4(n) describe the direct (coherent) 
and diffused (incoherent) parts of the light field in a random medium, respectively. The 
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direction of the vector ño coincides with the direction of an incident beam. It follows from 
equations (17) and (18) that 


Èe (iio) + Fal) = -N ÊF, (iio) — NLFa(i) +N | dy fA’, A) Fa’) S RA) 

+N f (ño, n) Fe (ño) f* (io, n) + BU) (19) 
where 

iy dF ee 
F(n) = —=(1-V)F. 
dn 
Thus, one can obtain 

F. (iio) = —NLF. (io) (20) 
for the coherent component and 
Fy(n) = —NLFy(i) +N / dy fR’, n) Fan’) f, n) + Bo(n) + Bn) (21) 
for the diffused light field. The tensor 


Bo(n) = Nf (io, n) Fe (fio) f* (Mo, ñ) (22) 


describes the single light scattering of an incident light beam. It can be easily obtained after 
calculation of F,(ño) from equation (20). 

It is well known that the coherent field A propagating in the direction ño satisfies the 
following formula [7]: 


E. (iio) = iANf (iio, fio) Ec (iio) (23) 


where A is the wavelength. 
From equation (23) and the definition of the coherent component of the light beam tensor 


Fo =) BES (24) 
one obtains 
Fe= X ESBS + EEY =iaN Y (JEES — ES f* ES) (25) 
or 
Fo =iAN(f Fo — Ff”). (26) 
It follows from equations (20) and (26) for the operator L that 
Ê Iy) = ià (flv) — (WI f+) (27) 


where we have used Dirac notation. 
Finally, we obtain the covariant radiative transfer equation (CRTE) for the diffused light 
(see equations (21) and (22)) 


G- OFA, A) = iAN [f G, DRED — FED fA, A] 
iN f E A aA 
+Nf (io, 8) Fo, Ro) f Go, R) + BO) (28) 


where F is the radius vector of the observation point. 
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The value of F.(7, no) in equation (28) is determined from (26), which can be written in 
the following form: 


(i - V) Fe = iàN ( f (fio, io) Fe — Fe f* io, fio)) (29) 
or 
F.(F, ño) = exp{—NhL} Fo(F, ño) (30) 


where Fo(/o) is the beam tensor for the incident light and A is the length of a beam path in a 
medium, which we consider to be uniform. Equations (27) and (30) provide the generalization 
of the Bouger law. 

The boundary conditions for equation (28) state that there is no diffused light arriving at 
a scattering convex medium from outside, 


Fafo, i) = 0 (at 7-1 <0) (31) 
where / is the unit vector normal to the boundary in the outward direction at the point with the 
radius vector Fo. 


Note that the general (without introducing separate coherent and diffused components) 
CRTE can be obtained from equations (17) and (27): 


G- V)F@, a) = iAN[ f(a, nF, n) — F(r,n) ft, n)] 


+N T E a a + BM (32) 
with the boundary condition 
F(ro,n) = Fo(Fo, n) (at 7-1 <0) 


where Fo is the radius vector of the point at the boundary of a scattering medium. 
Equation (32) can be transformed into the standard vector radiative transfer equation, 
assuming that 7 ||Z. It follows from equation (32) that in this case 


Gi. Oo, ñ) =iAN [AG * PE. ñ) — pF, AG, 1] 
+N f d AG DoE. RAGDE DO (33) 


where p is the density matrix (see equation (4)) and the matrix A describes the transformation 
of the electric field E = E NA + Eb due to light interaction with an elementary volume of 
a scattering medium in the coordinate system attached to a light-scattering layer. The matrix 
D(n) can be obtained from the source term B (7) in equation (32), assuming that 7 ||Z. It follows 
after substitution of equation (5) into equation (33) and simple algebraic calculations that 


Gi VIG, ñ) = —e TF, À) — enO, A) — enU C, A) — enV GF, ñ) 


on f 69% lol, AN) top OF, n) to,3U (7,7) + o4V(F, n')] + D 


van 
= 


n-V)OG, i) = —en I F, ñ) — en OG, A) — enU CF, A) — enV (F, ii) 


+N / dg [on 1 F,R) + on OG, n) + 023U (F, n) + or4V (F, 1')] + Do 


ye te Sy SA. SAn D IA 4.3 (34) 
-W)U(r,n) = —€3,1 (F, n) — €320 (r, n) — £€33U (r, Nn) — €34V (r, n) 


van 
= 


+N f dgy [031I (F, R’) +0320, n') +033U F, n’) + 034V C, n’)] + Dg 


n- V)V E, A) = —eqT(F, A) — eQC, A) — enU F, A) — enV E, n) 


van 
= 


+N $ dg [oul F, R’) +0 QC, n') + oUF, N') + oV C, n')]+ Dy 
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where expressions for elements oj; , €;; are presented in the appendix and D; = D,\+D2, D> = 
Dy; — Dz, D3 = Dj2+ D2, Da = i(Di2 — D21). Equations (34) can be rewritten in the vector 
form 


G- VSE, ñ) = -850P +N f a2 656. i’) + DCE, ñ) (35) 


where DF , n) is the vector with elements Di, Do, Ds, D, and matrices ê, ô have elements 
€ij, oij (i, j = 1,2,3,4), respectively. Equation (35) coincides with the standard vector 
radiative transfer equation [1]. It follows from this equation, within the framework of the 
scalar approximation, that 


(i -VIE ñ) =—e1G, iy +N f aro DI G, ii’) + BG) (36) 


where € = £11, 0 = 011. This is the well known scalar radiative transfer equation [2]. 


4. Conclusion 


In conclusion, we underline that the derived equations (28) and (32) can be used to study 
the polarization characteristics of a light field in different complex random media, including 
gyrotropic and anisotropic ones, within the framework of the coordinate-free approach. They 
also could be of help in derivations of simple and compact tensor expressions for isotropic 
media. 

The description of the interaction of light with matter (reflection, refraction and scattering 
of waves) in terms of the light beam tensor provides the possibility of considering the 
corresponding processes by covariant tensor methods, which in many cases simplifies the 
analytical derivations considerably. This was shown by Fedorov [4, 5] for the case of uniform 
anisotropic and gyrotropic media, where the selection of a coordinate system either related to 
the geometry of the problem (e.g. a plane-parallel layer) or to the symmetry of the dielectric 
tensor leads to extremely cumbersome results. It is expected that the same progress can 
be achieved for random complex media. Equations (32), (33) and (35) differ due either to 
the different coordinate systems or to the different description of the light beam. However, 
they provide the same information in the end, namely the intensity of the light, the degree 
of polarization and the characteristics of the polarization ellipse of reflected, transmitted and 
internal light fields. The choice of which equation to solve is largely dependent on the problem 
at hand. 

For instance, radiative transfer in isotropic random media (e.g. fogs and clouds) can be 
studied with equation (35). Equation (32) simplifies radiative transport studies in complex 
media (e.g. asymmetric and anisotropic ones). 

To solve equation (32) one needs to know the 3 x 3 scattering tensor f (n, n’) of a random 
medium which depends on the size of the scatterers, their shape, the dielectric tensors of the 
particles, etc. This tensor can vary with the location inside a medium. Many particles can 
contribute to the value of f(n’,n). Thus, the value of f(n, n’) is the average value of the 
scattering tensor for the ensemble of scatterers. 
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Appendix. Extinction and scattering matrices 


Elements of extinction é and scattering G matrices of the vector radiative transfer equation (35) 
are related to elements of matrix A in equation (33) with the following formulae: 


E11 = €22 = £33 = £44 = ÀN (A1, (0) + A, (0)) 
£12 = £21 = ÀN (A1, (0) — A5,(0)) 

£13 = €31 = ÀN (Ah (0) + Aj, (0)) 

E14 = £41 = ÀN (— Ah (0) + A}, (0)) 


£23 = —€32 = ÀN (Ah (0) — Aj, (0)) 

£04 = —£42 = —À N (Af (0) + A}, (0)) 

£34 = —£43 = ÀN (A3 (0) — Aj, (0)) 

o = 3 (lAn? + A21? + A121? + |4211") 


1 2 2 2 2 
on = į (lAn — Anl? + [Aail? — |Az2l’) 
Part / / / / H A Wt H 
013 = Ay Ayn + Ay Any + An An + A1412 
=< / ” n / n / / ” 
O14 = -A23 Ag, t Ag Ags + AY Ay AA 


oz = 4(-l42[ + |An + |Ai2|? — |Aail’) 


2 2 2 2 
on = § (|Az2|* + [Anl — JAn? — |Aail*) 


073 = Al Ah — A5 A) — AQ Ay + Ai Ain 
024 = —Aj Aly + Al Ain — AZ Ady + Ad) Ad, 
031 = AA}, + AAI, + AQ AT + ADA 
032 = —AH Al + ASAT) + AA) — An Aly 


= á id H H £ / / f 
033 = An Aj, + An A12 + An Aq + Ad) Ayn 


034 = A AY, — Ad Ai, — A Ah + A A12 
04, = -AA + A Ay — AAT, + AD Ai, 
oy = Ay Aly — An Ain — AAT, + ANAL 
043 = ~A Al, + ADA, — A Aly + AD Ain 
O44 = [ Ai, + AD AT — Ad Ain — Adi 12] 


where A’ j = Re(Aij), Aj, = Im(A;;). The symbol A;; (0) denotes the value of the amplitude 
matrix element in the forward scattering direction. 
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